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Abstract: The probability distribution D{e) in the energy loss incurred by incoming 
and outgoing hard quarks in a QCD medium is computed numerically from the 
BDMPS gluon spectrum. It is shown to follow an empirical log-normal behavior 
which allows us to give the quenching weight a simple analytic parameterization. 
The dependence of our results under the infrared and ultraviolet sensitivity of the 
gluon spectrum is investigated as well. Finally, as an illustration, we discuss and 
compare estimates for the quenching of hadron spectra in nuclear matter and in a 
quark-gluon plasma to HERA and RHIC preliminary data. 
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1. Introduction 

The energy loss experienced by a fast parton may serve as a measure of the density 
of color charges of the QCD medium it travels through [1, 2]. It may be huge in 
a dense medium such as the quark-gluon plasma (QGP), expected to be formed in 
the early stage of ultrarelativistic heavy ion collisions. One of the most direct and 
observable consequences would be the quenching of high particle production [3, 4]. 
This is reported experimentally. The first PHENIX data showed a large suppression 
in central Au-Au collisions as compared to proton-proton measurements [5], later 
confirmed on a wider p± range by the PHENIX [6], PHOBOS [7], and STAR [8] col- 
laboration. The quenching factor of Q{p±) ~ 5 observed at large p± finds therefore 
a natural explanation in this context. This makes possible tomographic investiga- 
tions [9], i.e., the study of the (color) structure of hot QCD matter by measuring the 
energy attenuation of hard probes. 
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Although the understanding of medium-induced parton energy loss has been 
extensively developed over the last few years [1, 2, 10, 11, 12, 13, 14, 15], less is known 
about how to relate this mechanism to observable quantities. A step in that direction 
has however been taken recently by Baier, Dokshitzer, Mueller, and Schiff (BDMS) 
in Ref. [16] in which they connect the measured quenching factor Q{p±) to the 
induced gluon spectrum dl/duj radiated by the leading parton. Perhaps even more 
importantly, the authors emphasize that the standard modeling of the quenching 
- determined by the mean energy loss - proves inadequate. Rather, the knowledge 
of the full probability distribution D{e) in the energy loss is actually required. On 
general grounds, the cross sections will be modified in the medium as 

^medium D{(l) ® ^vacuum (11) 



However, unless strong assumptions are made as for dl/dcu, D{e) cannot be calculated 
analytically in a simple way. It is therefore the aim of this paper to present and 
discuss a numerical computation of the distribution D{e) from the medium-induced 
gluon spectrum derived by Baier, Dokshitzer, Mueller, Peigne, and Schiff (BDMPS). 
Using these results, we give some estimates for the quenching ("diagnosis") to be 
compared to experimental data in electron-nucleus and nucleus-nucleus collisions. 

The outhne of the paper is as follows. In Section 2, we detail the numerical 
procedure to compute D{e) from its integral representation given in [16]. The dis- 
tribution is simply related to the gluon multiplicity N{ijj) radiated by the leading 
parton, calculated in Sect. 2.2. Results are given in Section 3. The probabihty dis- 
tributions for both an incoming and an outgoing quark are first presented (Sect. 3.1) 
and then contrasted with analytical approximations (Sect. 3.2). A simple analytical 
parameterization of the quenching weight D{e) ends the Section. We examine in 
Section 4 the sensitivity of our results beyond the approximations made in Section 3 
for the infrared and ultraviolet behavior of the BDMPS medium-induced gluon spec- 
trum. We apply our results in Section 5 where the quenching of hadron spectra is 
determined in cold and hot QCD matter and compared quahtatively to HERMES 
and PHENIX prehminary data. Finally, Section 6 summarizes the main results of 
the present study. 



2. Computation of the energy loss distribution D{e) 
2.1 Integral representation 

The multiple soft collisions undergone by a hard parton traveling through a medium 
induce gluon emission. Consequently, these radiated gluons take away an energy e 
from the leading particle with a probability distribution (or quenching weight) D{e). 
Let us briefly recall here its expression found by BDMS [16]. 
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As long as interference effects between radiated gluons (suppressed by as) can 
be neglected, we may assume that the gluon emissions from the leading parton are 
independent^. This allows for a Poisson formulation of D{e), which reads [16] 



^w-s:;^ n/ 

n=0 i=l '' 



dl{uji 
duj 



1=1 



■ exp 



duj 



■ (2.1) 



Here, dl/dui represents the medium-induced gluon spectrum and n the number of 
radiated gluons by the hard parton. Note that secondary gluon emissions, neglected 
in the soft limit, are not taken into account in (2.1). 

Using the Mellin representation of the delta function, the series (2.1) is resumed 
to finally obtain an integral representation 



D{e) 



dv ~ , ^ , 
2m ^ ' 



(2.2) 



with the integration contour C chosen (here) to be the imaginary axis. The Laplace 
transform D{'i') is simply related to the induced gluon spectrum through 



D{v) — exp 



-V 



f 

Jo 



duj e"^'" N iuj) 



(2.3) 



where the integrated gluon multiplicity A^(ci;) is defined as the number of gluons with 
an energy larger than a;, i.e., 

, dI{J) 



N{uj) 



duj 



duj' 



Taking u = ib in (2.2), the distribution D{e) thus becomes 

db 



D{e) 



with 



/ 

m 
m 



TT 



exp(-6/,(6)) cos(6(e-/e(6))) 



dcu cos (6a;) N{u!), 



(2.4) 



(2.5) 



+ 0O 



(2.6) 



diu sin (6 a;) A^(a;) 



Given a medium-induced gluon spectrum dl/diu, the integrated gluon multiplicity 

A^(a;) can be determined exactly. This then allows for the numerical computation of 
the integrals (2.6) and, subsequently, of the quenching weight D{e) through Eq. (2.5). 
The multiplicity of soft radiated gluons computed from the BDMPS spectrum for 
outgoing and incoming quarks is now discussed. 

"'^ Strictly speaking, one has to assume moreover that the energy loss e remains much smaller than 
the parton energy e <C (soft limit), which is precisely the case of interest in this section. We shall 
come back to this point in Section 4 where a finite quark energy will be explicitly considered. 
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2.2 Gluon multiplicities 



Outgoing quarks 

Let us start with the BDMPS gluon spectrum radiated by an outgoing quark 
with energy E traversing a medium of length L. In the soft gluon approximation 
(u; < £;), it reads [4, 12] 

= -InlcosUl + i)u]\ (2.7) 



in [ cosh u — sm u\ ; = \/ — ' 



2uj V 2 a; tt 

where Cp — 4/3 is the Casimir operator in the fundamental representation and 
as = g^/4:7r ~ 1/2 the strong coupling constant . In the soft limit, the spectrum (2.7) 
is thus characterized by only one energy scale 

uJc=^qL^ (2.8) 

where the so-called gluon transport coefficient q measures the "scattering power" of 
the medium. Relating it to the gluon density of the medium, BDMPS give pcrtur- 
bative estimates for the transport coefficient. While it is shown to be as small as 
q ~ 0.25 GeV/fm^ in cold nuclear matter, a much larger g ~ 5 GeV/fm^ is expected 
in a hot {T — 250 MeV) quark-gluon plasma [11]. 
Writing N{u) as 

N{uj)^ / duj' ; ' - / dw' ; ' (2.9) 

and Taylor expanding the gluon spectrum at^ a; = a;/a;c <^ 1, the gluon multiplicity 
reads [16] 

N{u<^Uc) - a |^y| + ln21ncl;- 1.44136 | (2.10) 

in the infrared (IR) domain, while it drops much faster 

all 



in the ultraviolet (UV) region. It is however necessary to rely on a numerical de- 
termination of N{ijj) in between these two regimes. We decompose therefore the 
integrals (2.6) as 







/"OO 

+ / 


+ y 


du! COS ipuj) N{u}) 


(2.12) 




. •/() 











^We shall keep this notation for all variables normalized to the energy scale ujc- 
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where the first two terms on the r.h.s. of (2.12) are expressed analytically in terms of 
Fresnel integrals and hyper geometric functions'^ while the third term is determined 
numerically. The cutoff ui^ (a;+) under (above) which the exact gluon multiplicity 
may be replaced by its analytical approximation (2.10) (respectively, (2.11)) is taken 
to be uj- — 0.01 = 2u!c)- 

The gluon multiplicity radiated by an outgoing quark N{u!) is plotted in Figure 1 
as a function of Cj {upper solidlme), together with its analytic approximation (2.10) at 
small energy (dashed). Anticipating the discussion in Section 3.2, the approximation 



N{uj) = a\- (2.13) 
V uj 

based on the small Co behavior of (2.10) is shown as a dash-dotted line. While the 
expression (2.10) is shown to reproduce the exact result for gluon energies up to 
CO ~ 0.3 a;c, Eq. (2.13) strongly overestimates N{u!) over the whole range of interest. 

Incoming quarks 

A similar procedure for the case of an incoming quark can be carried out. The 
medium- induced gluon spectrum now reads [12] 

sin [ (1 + i)u 



dliuj) a , 
— T^-^ = — m 
duj uj 



= ^ln 

2uj 



{l + i)u 
cosh^ u — cos^ u 



(2.14) 



2^2 

As previously, the limiting energy behavior can be extracted analytically. We found 




N(u;^u;c)^a \ ^/ - + In 2 IncD - - (Inu;) - 1.32099 | . (2.15) 
and 

at small and high energies respectively. The integrals (2.6) are computed as for the 
outgoing quark case. 

The energy dependence of the gluon multiplicity emitted by the incoming quark 
is also shown in Figure 1 (lower solid line) together with the analytical expression 
(2.15) as a dashed line. The number of emitted gluons remains much smaller than 
what is observed for the outgoing quarks (note the factor 15 difference between 
(2.11) and (2.16)). Indeed, for a hard quark produced inthe medium, a gluon can be 
emitted shortly after the hard process and prior the first scattering of the hard quark 
in the medium [12, 17]. This increases the gluon multiplicity (hence, the energy loss) 
radiated by outgoing quarks prod uced in the medium. 

^Even though the gluon multipHcity in the BDMPS formaUsm is sensitive in the IR, all integrals 
here are infrared safe quantities. 
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co/co 

c 

Figure 1: Integrated gluon multiplicity N(uj) eomputed from the BDMPS spectrum re- 
spectively for outgoing (upper solid) and incoming (lower solid) quarks. Its analytic IR 
behavior (resp. Eq. (2.10) and (2.15)) is shown in dashed lines while the dash-dotted line 
represents the small u approximation (2.13), N{uj) oc 

3. Numerical results 

3.1 Incoming vs. outgoing quarks 

Following the procedure detailed in the former Section, the probability -D(e) that an 
outgoing or incoming quark loses an energy e while going through the medium is 
determined numerically. 

The normalized distribution 

D{e = e/uj,)=uj,D{e) (3.1) 

is represented in Figure 2 for the outgoing (solid line) and the incoming (dashed) 
quark case. It exhibits a strong peak for energy loss e much smaller than the typical 
scale ujc- Moreover, the long energy tail of the distribution makes the mean energy 
loss (e) 

(e) = / deeD(e) (3.2) 
Jo 

well larger than its peak value, as emphasized in [16] . 

Furthermore, Figure 2 shows that the energy loss suffered by an incoming quark 
proves much smaller than what happens to a quark produced in the medium. This 
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Figure 2: Distribution in the energy loss D(e) for outgoing (solid) and incoming (dashed) 
quarks computed from the BDMPS spectrum in the soft hmit. 

was already apparent from the gluon multiplicities (Figure 1 and end of Sect. 2.2): 
the lower the number of emitted gluons for all cu, the smaller the energy loss e. 

3.2 Comparison with the BDMS estimate 

The numerical results we obtained previously are compared with analytical estimates 
for the quenching weight. Let us begin with the illustrative guess based on the small 
energy behavior of the gluon multiplicity (2.13). Using Eq. (2.13) in (2.3), the inverse 
Laplace transform reads 

D(i') = exp (^—a^/27^PUJ^ . (3.3) 
which leads, via (2.2), to [16] 

The expression (3.4) is plotted in Figure 3 (dashed) as a function of the energy 
loss e, together with the full calculation for outgoing quarks (solid). Whereas the 
location of the peak of the distribution is roughly similar (e ~ 0.1), it exhibits a much 
larger energy tail than what we found numerically^. Consequently, the probability 
for a small euerg}' loss is soniowliat reduced as compared to our full result. 
^Remark for instance the infinite mean energy loss (e) in Eq. (3.4). 
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^ 3.5 




Figure 3: Comparison of the distribution -D(e) computed numerically for outgoing quarks 
with the BDMS analytic estimate (3.4) from the small u) behavior (2.13). Also shown in 
dotted line the analytic formula for D{e) extracted from the approximated gluon multiplic- 
ity (2.10). 



The poor agreement between the illustrative guess (3.4) and our exact prob- 
ability distribution could have been anticipated from Figure 1 where the small uo 
approximation for N{uj) {dash-dotted) was shown to fail on the whole energy range. 
On the contrary, as stressed in the previous section, one might guess that the gluon 
multiplicity (2.10) (respectively, (2.15)) emitted by an outgoing (respectively, incom- 
ing) quark gives more satisfactory results. The inverse Laplace transform D{v) now 
reads [16] 

D{v) ~ exp [-a y^huJu^ - In 2 \^{vu^ - 1.84146)] , (3.5) 
for an outgoing quark case, whereas we found 



i){y) — exp 



-a 



^yJ^huJu^ - 0.981755 \\i{yu^ - ^ \x^{vu:^ - 2.21561^ 



(3.6) 



for incoming quarks. Therefore, the quenching weights D°"*(e) and D^'^ie) can be 
determined analytically through the inverse Laplace transforms of (3.5) and (3.6), 
respectively. 

Although analytic, the expressions obtained are lengthy and hardly transparent 
(sum of hypergeometrical functions), and are thus not reproduced here. Rather, we 
display in Figure 3 the excellent agreement between the analytic formula for D°"*(e) 
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(dotted) together with the result computed from the exact gluon spectrum. We may 
notice in particular that the agreement remains perfect up to energy loss e ^ 0.3 c^c, 
above which the small energy approximations (2.10) and (2.15) for A^(a;) start to fail 
(See Section 2.2 and Figure 1). 

3.3 Analytical parameterization 

As mentioned in the introduction, the probability distribution bridges the gap be- 
tween the theory of medium-induced parton energy loss on the one hand and the 
observable consequences on the other hand. In particular, its knowledge is required 
to model the quenching of hadron spectra in nuclear collisions. This was the main 
motivation for our present study. 

However, we have seen that D{e) cannot be solved analytically, and thus neither 
can the medium cross section (1.1). Even though we have just stressed that an 
analytic expression is shown to mimic almost perfectly the numerical results, its 
complicated expression makes it useless in practical terms. For the sake of simplicity, 
we shall give instead in this section an empirical analytical expression for D{e). 

The energy dependence of the quenching weight £)(e) follows with a great accu- 
racy a log-normal distribution. 



1 



exp 



(loge- fiY 
2a2 



(3.7) 



characterized by two parameters, /i and a. The very nice agreement between the 
log-normal parameterization and the full result is illustrated in Figure 4 for outgoing 
(left) and incoming (right) quarks, using the parameters given in Table 1^. 





out 


in 




-1.5 


-2.55 


a 


0.73 


0.57 




0.3 oJc 


0.1 LOc 



Table 1: {^,<t) parameters of the analytic approximation (3.7) to the distribution D{e) 
for both outgoing {left) and incoming {right) quarks. 

The simple analytic expression (3.7) can therefore easily be used to compute the 
medium cross section (1.1) and hence the quenching of hadron yields. 

^Although wc give this two parameter set, a small but clear anticorrelation exists between and 

a. 
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To give the reader another feehng for the agreement between the parameteriza- 
tion (3.7) and the exact result, the mean energy loss experienced by the fast parton 
in the medium, 

(e) = exp (^^^ , (3.8) 

is computed in Tabic 1. The values tend to be pretty close to the result found 
analytically by EDMS [12] 

(e) = ujc = ^ ojc, (outgoing quark) 

(e) = — - — ujc^ -ujc- (mcommg quark) 

D y 

EDMS [16] emphasized that the standard modeling of the quenching using the 
mean energy loss 

^medium _ ^(^^)) ^ ^vacuum (g^^^Q) 

instead of (1.1) leads to a strong bias — when cr^'^'^""™ is a steeply falling function — 
whose strength is given by the higher moments of ^(e) [16]. In particular, the larger 
the asymmetry of the quenching weight (characterized by its skewness parameter 
7 = (e^)/(e^)^/^), the stronger the bias effect. Indeed, we remark for instance in 
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Figure 4 that the mean (e) proves larger than the most probable value of the distri- 
bution. In addition to that, let us note that this bias will become more pronounced 
when considering outgoing than incoming quarks, from the larger skewness of the 
distribution (7 ^ 3.1 and 2.1, respectively). 

Although certainly useful, one may nevertheless wonder about the origin of the 
log- normal dependence of the probabihty distribution -D(e). The reason for this 
is not clear. Even though it could be accidental, we argue it may rather come 
from the effect of the combinatorics in the Poisson approximation (2.1). Let us be 
more explicit. Because the BDMPS medium- induced spectrum dl/duj dramatically 
drops with the gluon energy uj, the energy loss e will be carried away by a large^ 
number of soft radiated gluons. It could be the huge product of radiation probabilities 
in the Poisson approximation (2.1) that is responsible for such a behavior. We 
therefore conjecture that a log-normal quenching weight may emerge for any fast 
falling medium-induced gluon spectrum. This will be further discussed in the next 
section where the series (2.1) is partially summed. 

4. Approximations 
4.1 Bethe-Heitler regime 

The lifetime t of the medium induced gluons emitted by the hard quark is given 



where the transverse momentum k'j_ of the gluon and its mean free path Xg are 
related through the transport coefficient q = fJ^'^/Xg [12]. As long as the lifetime t 
remains much smaller than the typical distance between two scattering centers in 
the medium, the radiation spectrum will be proportional to the gluon emission on 
one single radiator, that is 



where L/Xgis the number of collisions encountered by the leading quark with a mean 
free path Xg. This is the Beithe-Heitler (BH) regime that occurs in QCD for small 
lifetime gluons. 

On the contrary, a gluon with a long lifetime t (as compared to its mean free 
path Xg) will only see a group of scattering centers as a whole. In this regime, the 
induced gluon spectrum (4.2) is suppressed by the number of scattering centers Ncoh 

^Strictly speaking, the series has to be entirely resummed because of the vanishing term 
exp(-Ar(0)) in Eq. (2.1). 



by [2] 




(4.1) 




(4.2) 
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that act coherently [2], 



d7(u;)™_ L fdl{uj)\ 1 



duj Xn \ duj J Ncoh 



= - X 1 ) X (4-3) 

(1) 



which is the Landau- Pomeranchuk-Migdal (LPM) gluon spectrum (2.7) and (2.14) 
computed by BDMPS. Hence, this coherent regime will set in for gluon energies 
greater than 

l^min = k^Xg^ jl'^ Xg. (4.4) 

A rough estimate Umin ~ 300 MeV for a hot QCD medium is given by EDMS in [16]. 
At small gluon energy, the induced gluon spectrum is negligible as compared that of 
the extrapolation of the LPM spectrum (2.14), 

dl(uj) d7(a;)™ 
cu — 7^ = constant <^ uj — at a; < uJmin- (4.5) 
acu duj 

As already pointed out at the end of the previous section, the quark energy loss 
actually originates from the emission of a large number of soft gluons, with en- 
ergy a; ^ e ~ (9 (cJc), because of the divergence of the BDMPS gluon spectrum in 
the infrared sector. Hence, we anticipate a large sensitivity in our results to the 
Bethe-Heitler regime (a; < uJmin)- To estimate this sensitivity, we have repeated the 
computation of the quenching weight using a somewhat arbitrary truncated induced 
gluon spectrum, 

dl{u < u^^) ^ ^ 
dijj 

where the IR cutoff uJmm is expressed as a function of the scale ujc- 

As a consequence, the number of gluons radiated by the hard quark A^(0) = 
N{uJmin) becomes finite, unlike in Section 2 where it was shown to be slowly divergent 
in the IR. Using (4.4) in (2.13), it is given by the number of collisions L/Xq (or 
opacity) in the medium [16]. Therefore, the probability po for having no interaction. 



Po = exp 



du 



exp(-L/A,) (4.7) 



that appears in Eq. (2.1) no longer vanishes in a finite length medium. In the general 
case, the probability distribution will read [18, 19] 

D{e)=po6{e) + d{e). (4.8) 

Following Wiedemann [18], we shall subtract the discrete contribution in (4.8) before 
performing the numerical inverse Laplace transform, 

d{e) = j^^.{b{v)-p,)e- (4.9) 
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to eventually extract the continuous part d{e) of the probability distribution. 

The result for d{e) assuming cUmin — is shown in Figure 5 (dashed) to- 

gether with the previously computed D{e) without any IR cutoff (solid). Taking 
g ~ 5 GeV/fm^ and L = 5 fm in (2.8), this would correspond to cUmin — 600 MeV. 
First, Figure 5 clearly indicates that the shape of the distribution looks pretty similar 
to what has been obtained before, although shifted to smaller energy loss. This does 
not come as a surprise as fewer gluons are radiated. Another remarkable feature is 
the structure observed at e <^ lOc- This actually originates from the emission of a 
very small number of gluons n in the Poisson series (2.1). To go a bit further, we plot 
in Figure 6 the quenching weight as a function of e in units of the cutoff cUmin- The 
distribution (solid line) is identically equal to zero up to e = comin when the channel 
for the one gluon emission opens, followed by a strong decrease coming from the 
dropping gluon spectrum (2.7). Angular points may clearly be seen at t/umin = 2 
and 3 which correspond to the opening of the two and three gluon radiation channels, 
respectively. 

Because each term is finite in the Poisson expression (recall, po > 0), it becomes 
possible to sum all the individual contributions to reconstruct the quenching weight, 
as 

m— 1 

D(e < muj^,n) = PoS(e) + ^^"H^), (4-10) 

n=l 
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Figure 6: Distribution d(e) for outgoing quarks computed from the BDMPS spectrum with 
an IR cutoff comin = t<Jc/100 (solid) together with the probabihty distributions u}cd^'^\e) to 
have i = 1 (dashed), 2 (dotted), 3 (dash- dotted) gluons emitted (see text for details). 



where the S"'' represent the probabihty for having exactly n gluons emitted and 
where the terms n > m ^ ^/oJmin do not contribute. Although the resummation 
procedure offered by BDMS allows one to compute elegantly all these terms, we 
nevertheless compute the first three terms of Eq. (4.10), which can be written as 

<(<'>(£) =Po X ^ X e(e-w„i„), (4.11) 
= I X r-'- d.. '-1^ X e(. - 2.„.„), (4.12) 

and, 

(4.13) 

for the emission of 1, 2, and 3 gluons. These are displayed in Figure 6 as dashed, 
dotted, dash-dotted lines, respectively. Their sum is shown to reproduce exactly 
the full result in solid line. At larger energy loss e ^ ujmin, the sum over the large 
number m ~ e/cumin ^ 1 of gluons makes the distribution d(e) a smooth function of 
the energy e. We then recover the log-normal behavior previously discussed. 
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4.2 Beyond the soft gluon approximation 

The calculations performed so far have been computed in the soft gluon approxima- 
tion, i.e., assuming the gluon energy u to remain small with respect to the leading 
quark energy E. This approximation is, however, rarely justified in practice. BDMS 
give in Ref. [12] the full expression of the induced spectrum radiated by a quark. 
Neglecting O {{u/EY) terms, it can be written as 

where the additional factor actually comes from the quark-gluon DGLAP splitting 
function. 

The probability distribution D{e,E) is represented in Figure 7 for various quark 
energies E/uJc- The effect of the O (uj/E) corrections in the gluon spectrum (4.14) 
is to reduce hard gluon emission, and hence the high-energy tail of the distribution. 
On the contrary, the small e behavior of the quenching weight remains unchanged, 
with the exception of the absolute magnitude which follows from the normalization 
constraint. In particular, the location of the peak does not exhibit a strong quark 
energy dependence. Perhaps more interesting is the following observation. Although 
Eq. (4.14) ensures that a single gluon cannot carry more energy than available {cu < 
E), the use of the spectrum (4.14) does not a priori guarantee the quark total energy 
loss 

e = ^ c^i (4.15) 

i 

in the Poisson expression (2.1) to be bounded. Indeed, we do observe in Figure 7 a 
small but significant contribution of the probability distribution in this kinematically 
forbidden region, D{e > E) ^ 0, when the quark energy E is small enough. This 
signals the breakdown of the eikonal approximation on which the BDMPS framework 
rely. It is indeed no longer justified to consider multiple successive and independent 
quark-nucleon scatterings when the quark energy is smaller than, say, half the energy 
scale ujc- 

The need for a simple analytic parameterization of the probability distribution to 
compute the quenching of hadronic spectra has been stressed in Section 3.3. Hence, 
we shall now extend this parameterization for any quark energy. Noticing that the 
distributions D{e,E) still follow log- normal distributions, the quark energy depen- 
dence will enter through the parameters /i and a. The quenching weight thus reads 

(log6-^(^))' 



D(e, E) = \ exp 



2a{Ef 

where the ^Ji{E) and (7{E) are given by the empirical laws 

^l{E) = -1.5 + 0.81 X (exp (-0.2/^) - l) , 



(4.16) 



(4.17) 

a{E) = 0.72 + 0.33 x (exp {-Q.2/E) - l) . 
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Figure 7: Outgoing quark energy dependence of the distribution D[e, E) computed from 
the BDMPS spectrum Eq. (4.14). Shown in sohd hne is the previous result obtained in the 
soft Hmit. 



Using the parameterization (4.16), the mean energy loss (e) is plotted as a func- 
tion of the quark energy [dashed) in Figure 8. In Ref. [4], BDMPS discussed in a 
schematic way the characteristics of the E dependence of the mean energy loss (e) . 
They found that [4] 

^')/^^^y 1 if^>., (4-18) 



These features are quahtatively reproduced here. At high energy E ^ Uc, (e) is 
independent of E and rather close to its asymptotic value (e) ~ 0.3 (shown in solid 
line), whereas a stronger dependence is seen at low energy. E ~ oOc- In particular, it 
is interesting to note that the dependence of (e) will not set in until approximately 
^ 3 o^c- Our result is also reminiscent of what has recently been found by Gyulassy, 
Levai, and Vitev in the computation of the mean multiplicity of radiated gluons [20], 
that strongly increases up to a quark energy E ^ ujc ^ 15 GeV in their calculation. 
Finally, we show in Figure 8 {dotted) the expected energy loss (that we define as the 
most probable value of the quenching weight) that is well smaller than the mean (e) 
and independent of the quark energy. 



- 16 - 



W 0.4 



— mean 



• • • expected 



0.3 



0.2 



0.1 











1 



2 



3 



4 



5 



E/co. 



c 



Figure 8: Mean (e) (dashed) versus expected (dotted) energy loss as a function of the 
outgoing quark energy E. The sohd hne represents the mean energy loss determined in the 
soft limit. 

5. Applications 

The probability distributions for both incoming and outgoing quarks have been com- 
puted in the previous sections. To illustrate the use of these results, we now determine 
the quenching of hadron spectra in nuclear collisions and compare it to experimental 
preliminary data. 

5.1 Medium modification of fragmentation functions 

The HERMES collaboration at DESY recently reported on hadron yields measured 
in electron-nucleus collisions. They measured the production ratio 



in a "heavy" (N and Kr) over a light (D) nucleus for a given hadron species h. Here, 
u denotes the virtual photon energy in the lab frame, z the momentum fraction 
carried by the produced hadron, and where the multiplicity of produced electrons 
normalizes the hadron yield N^. 

The hadron multiplicity in (5.1) can be computed perturbatively to leading order 
(LO) in as- It is written in terms of (nuclear) parton densities qf{x,Q^,A) and 




(5.1) 
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fragmentation functions (FF) Dl!^{z,Q'^, A), 

j dxY,e)<lfM^A)a^*\x,v)D){z,Q\A) j (5.2) 

y dxY,e)q^{x,Q\A)a^*\x,v) 



1 dN\ 
N% dv dz 



and where the LO cross section is given by 

(5.3) 



X s 



The integral over Bjorken x — Q'^ji^Mv) appearing in (5.2) is given by the Q"^ 
acceptance of the HERMES experiment. To a first approximation, only the valence 
up quark will contribute to the hadron yield (5.2) when x is not too small. Hence, 
the ratio (5.1) will approximately be given by the ratio of the u ^ h fragmentation 
functions 

R\{z, v) ~ D^(;., Q\ A) I D^^{z, Q\ D). (5.4) 

Therefore, the nuclear dependence of the fragmentation functions might be revealed 
through the measure of B!^. We further note that the effects of nuclear shadowing 
in the parton densities qf{x, Q^, A) should remain small as they mainly cancel in the 
ratio N^/N^ (5.2). 

The multiple scattering of the produced quark in the nuclear medium may be 
responsible for the observed dependence of the ratio R\ with v and z [22]. Indeed, 
the energy loss of the hard quark will shift its energy from ~ i/ to E'^ ~ — e at 
the time of the hadronization. Therefore, this mechanism leads to a shift in z, 

- .• = ^ = -^. (5.5) 

u V — e 1 — e/v 

The effect of the energy loss mechanism on the nuclear fragmentation functions may 
be modeled according to [23, 22] 

zD'}{z,Q\A)= [ de D{€,u) z* D'}{z*,Q^). (5.6) 
Jo 

Using (5.6) in (5.2) and integrating over z, the u dependence of Rl^^v) in a 
krypton over a deuterium target has been computed. In the calculation, the energy 
scale ujc was determined by the transport coefficient q adjusted to the data and the 
length of matter covered by the hard quark proportional to the nuclear radius 



^Isospin corrections should be small as we compare nuclei with a similar Z/A ratio, and have 
thus been neglected. 
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Figure 9: Attenuation ratio R\^{v) (5.1) plotted as a function of u for charged hadrons 
(upper left), charged pions (upper right), positive and negative (lower) kaons. Calculations 
(solid) are compared to HERMES preliminary data (circles) taken from Ref. [21]. 



L = 3/4 R. The parton densities were given by GRV98 LO [24] while we made use 
of the Kretzer LO parameterization for the fragmentation functions [25]. 

The calculations for charged hadrons, pions, and kaons are compared with HER- 
MES preliminary data in Figure 9. The trend is reproduced well for all hadron 
species, although the calculation for the pions (7r"'" + 7r^) somehow underpredicts the 
effect. It is also interesting to notice that the K~ yield is more suppressed than the 
K'^ — as seen in the data — which arises from the stronger slope of the FF at large 
z. A much smaller difference between the tt"*" and 7r~ attenuation (average in the 
Figure) is also observed, unlike the present data which do not exhibit any isospin 



- 19 - 



dependence in this channel. Let us also mention that the z dependence of the ra- 
tio is fairly reproduced when z is not too large, the calculations predicting a much 
stronger suppression at large z than what is observed experimentally. However, the 
formula (5.6) we use may no longer be valid in this specific kinematic region where 
large higher twist corrections come into play [22]. 

The calculations were performed assuming the transport coefficient for nuclear 
matter to be g = 0.75 GeV/fm ^. This would correspond to a mean energy loss per 
unit length —dE/dz = {e)/L to be —dE/dz = 0.62 GeV/fm in a large (L f» 5 fm) 
nucleus. This result is close to what Wang and Wang determined in their analysis [22], 
although they had not considered the full probability distribution (but its mean) and 
neglect the u dependence of the mean energy loss (which is found to be relevant at 
low i/ ?a 8 — 10 GeV). Furthermore, let us remark that this estimate is reduced by a 
factor of three for incoming quarks, —dE/dz ~ 0.21 GeV/fm, in excellent agreement 
with our recent Drell-Yan data analysis [26] . 

It is not our aim to claim that energy loss is the only mechanism responsible for 
the trend observed in HERMES data. Indeed, many other effects such as nuclear 
absorption [27], gluon bremsstrahlung [28], or partial deconfinement [29] have been 
advanced to account for these measurements. Hence, the value extracted for q has 
to be seen as an upper limit. 

5.2 Quenching of high p±_ particles 

A huge parton energy loss is a clear signal for quark-gluon plasma (QGP) formation. 
Therefore, the production of high hadrons might reveal the existence, and more 
importantly the characteristics, of the produced medium in heavy ion collisions. 
Interpreting the depletion of the hadron yield observed at RHIC as coming from the 
energy loss experienced by hard quarks in an expanding QGP, we determine and 
compare the quenching factor 

= 7^^:^ (5.7) 

to PHENIX preliminary data on vr*^ production in central Au-Au collisions at ^/snn = 
200 GeV. The quenching factor (5.7) may be determined through the ratio^ [16] 

where the p_L differential vacuum cross sections can be computed perturbatively. Fol- 
lowing EDMS, we shall however adopt the fit proposed by the PHENIX collaboration 

J vacuum/ , A 



oc (1.71 + p^[GeV])-''-^* (5.9) 



*We neglect the longitudinal momentum of the produced quark (p ~ 
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for the sake of simplicity. To determine (5.7), we use our parameterization (4.16) for 
D{e, pj_), the energy scale cUc being given by the dynamical scaling law of Salgado 
and Wiedemann [19] 



dt'{t'-to)x i^jj . (5.10) 



where q{to) is the transport coefficient of the medium at the initial time to and a — 1 
characterizes the longitudinal expansion of the QGP. We take in the following an 
initial time to — 1 fm, a medium length L = 5 fm and the transport coefficient 
q{to) = 3.5 GeV/fm^ (which leads to (uc) ^ 10 GeV). 




Figure 10: left: The quenching factor Raa{p±) is computed from the BDMPS gluon 

spectrum in the soft limit with (dotted) and without (solid) IR gluon energy cutoff ojuiii-i. 
Effects of O {uj/p±) corrections to the soft gluon spectrum are shown in dashed Une. right 
: The influence on the modeling of the quenching in the p± dependence of the quenching 
factor is displayed. The solid line represents the calculation from the full probability 
distribution (5.8), while the quenching factor is computed assuming a shift of the mean 
[dashed) and the expected [dotted) energy loss in the cross sections (5.9). 



The quenching (5.8) is plotted as a function of p_L in the soft gluon approximation 
p_L ^ {ojc) in Figure 10 (solid). To give an idea of the uncertainties, we show on 
the left panel the effect of the O {uj/p±) corrections in the induced spectrum (4.14) 
(dashed), that slightly (~ 10%) reduce the depletion. The ratio Raa{p±) is also 
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calculated assuming a gluon energy cutoff uOmin ~ 100 McV {dotted). The right 
panel displays the inffuence of the modeling of the quenching. The suppression is 
computed by shifting the mean {dashed) and the expected {dotted) energy loss in the 
cross section (5.9) together with the result assuming the full distribution {solid). The 
effect of the bias is clear. The suppression is overestimated by more than a factor of 
two when one considers a mean energy loss rather than the full quenching weight. 
On the contrary, we note that shifting the cross sections by the peak value in the 
distribution underestimates the quenching. 

The quenching (5.8) is now compared to the PHENIX 7r° preliminary measure- 
ments [30] in Figure 11. For this illustration, the transport coefficient g(to) = 3.5 GeV/fm^ 
was chosen so as to reproduce high data^. Because of the strong correlation ex- 
isting between q and L for a given (a;c), this absolute value should however be taken 
with great care. Figure 11 clearly demonstrates that the trend of our estimate is 
opposite to what is observed experimentally, with the possible exception of the high- 
est data bins. It is not obvious to see where this discrepancy comes from. First, 
the extrapolation to low is somehow difficult, as our results shall be much more 
sensitive to the IR behavior of the BDMPS spectrum (2.7) where the LPM regime 
is no longer guaranteed. Secondly, we have already stressed that the eikonal approx- 
imation (2.1) breaks down for quark energies = O {{uJc)/2) ~ 5 GeV. Moreover, 
other mechanisms may compete and weaken the quenching in this low p± region, 
such as the detailed balance process where the leading quark picks up thermally 
equilibrated gluons in the QGP [31]. Let us finally mention that Vitev and Gyulassy 
recently suggested that the interplay of k± broadening, nuclear shadowing together 
with parton energy loss may account for the trend of the data [32] . 

It has been assumed so far that the hard quark propagates through the QGP 
with length L and subsequently hadronizes outside the medium. This picture should 
be true at high p± when the hadronization time is large enough because of Lorentz 
dilation. In the following, let us suppose that, in the pj_ range of interest here, the 
hard quark hadronizes inside the medium. The reader may of course worry about 
the relevance of the hadronization concept in a QGP. Thus, we rather imagine that 
hadronization might occur (and makes sense) in a cooling and more dilute {t ^ to) 
system such as a hot pion gas, although with a significant transport coefficient [33]. 
The length covered by the parton is then no longer given by the system size L, but 
its hadronization time th, 

th ^ K{z) X ^ (5.11) 
a 

where a ~ 1 GeV/fm is the string tension. Several models have been proposed in 
the literature to characterize the z dependence of the hadronization time [27, 28, 34]. 

^In this calculation, we made use of the parameterization (4.16) for D{e,p±_), i.e., with O {ui/p±) 
corrections and a vanishing gluon energy cutoff iVmin = 0. This corresponds to the dashed line in 
the left panel of Figure 10. 
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Here, we take K = 1/2 assuming a mean {z) = 1/2 cither in the approach by 
Brodsky and Mueller [34] or in the gluon brcmsstrahlung model [28]. In other words, 
the length ~ th over which the quark propagates in the medium actually depends on 
(and increases with) p±. We may write the mean (cUc) in this simplified picture as 



The formula (5.10) should be replaced by (5.12) which depends approximately lin- 
early with p± and thus so does the energy loss. As can be seen in Figure 11, this 
leads to a decrease of the quenching as a function of p± (as long as th{p±) — O (L) 
over which Eq. (5.12) ceases to be valid) in good agreement with the trend of the 
data assuming g(to) = 4.5 GeV/fm^ In this calculation, we have not considered 
possible final state interactions of the produced pion in the medium, which may 
somehow modify this picture. Provided these to remain rather small in a dilute sys- 
tem, a decrease followed by a subsequent increase of the quenching at larger p± would 
therefore be a signal of the transition between hadronization inside and outside of 
the hot medium. 

6. Conclusion 

Let us summarize what has been carried out here. 

The multiplicity of gluons N{u!) radiated by a hard quark in a QCD medium has 
first been calculated from the BDMPS gluon spectrum in the soft hmit {coc <^ E). 
The analytic IR behavior of A'"(a;) for outgoing quarks is recalled while we give similar 
expressions considering incoming quarks as well. These are shown to reproduce fairly 
well the exact multiplicity for gluon energies as large as the typical scale u ~ Uc- 

Subsequently, this allows for the computation of the probability distribution from 
the integral representation given in Ref. [16]. The quenching weight determined for 
both incoming and outgoing quarks is then compared to analytic estimates based on 
the small energy behavior of the gluon multiplicities. In particular, we have empha- 
sized that the analytic expression given by EDMS strongly overestimates the high 
energy tail of the energy loss distribution, while an analytic (although complicated) 
expression based on the multiphcity N{u; < (Vc) (Eq. (2.10)) reproduces fairly well 
the full result. Noticing that the quenching weight follows a log-normal distribution, 
we give D{e) a simple analytic parameterization. 

Going a step beyond in Section 4, the probability distribution is computed from 
the BDMPS spectrum truncated in the IR (uj > uj„iin) to ensure the LPM regime 

^°The quenching will essentially depend on the product K x g(to), both terms being poorly known. 
It is needless to repeat that this transport coefficient chosen to account for the data is thus a rough 
estimate. This prevents us from drawing any conclusion from its absolute value. 





<L (5.12) 
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Figure 11: Preliminary PHENIX data on vr suppression factor in Au-Au collisions at 
y/s = 200 GeV as a function ofp± (circles) compared to the theoretical calculations assum- 
ing hadronization inside (dotted) and outside (dashed) the hot QCD medium. Systematic 
errors are not shown. 



to be at work. It exhibits a "discrete" behavior at small energies (a few times the 
gluon energy cutoff) which corresponds to the emission of a small number of gluons, 
followed by a smoother "continuum" similar to what was found earlier. 

The 0(uj/E) corrections in the gluon spectrum are then exphcitly taken into 
account. This leads to a reduction of hard gluon emission which suppresses the tail of 
the probability distribution. To be complete, we give a (log-normal) parameterization 
for the quenching weight D(e, E) for finite quark energy, which can be of use for future 
tomographic studies. 

In conclusion, we illustrate the former results computing the quenching of hadron 
spectra in nuclear matter as well as in an expanding quark- gluon plasma. These esti- 
mates are respectively compared to HERMES e-A data and recent 7r° measurements 
by the PHENIX collaboration in Au-Au collisions at RHIC. While the HERMES 
data can quantitatively be understood as coming from the effect of quark energy 
loss, the trend observed in the PHENIX data is opposite to what one could naively 
expect. Finally, we suggest that this could be due to the fact that the "hard" quark 
produced at moderate p±_ actually hadronizes inside the medium. This has of course 
to be further investigated and is under current study. 
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